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We derive a computable analytical formula for the quantum fidelity between two arbitrary multimode Gaus¬ 
sian states which is simply expressed in terms of their first- and second-order statistical moments. We also show 
how such a formula can be written in terms of symplectic invariants and used to derive closed forms for a variety 
of basic quantities and tools, such as the Bures metric, the quantum Fisher information and various fidelity-based 
bounds. Our result can be used to extend the study of continuous-variable protocols, such as quantum telepor¬ 
tation and cloning, beyond the current one-mode or two-mode analyses, and paves the way to solve general 
problems in quantum metrology and quantum hypothesis testing with arbitrary multimode Gaussian resources. 
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I. INTRODUCTION 

The quantification of the similarity between two quantum 
states is a crucial issue in quantum information theory UIJ 
and, more generally, in the entire field of quantum physic sjlsll . 
Among the various notions, that of quantum fidelity [J-la] is 
perhaps the most well-known for its use as a quantifier of per¬ 
formance in a variety of quantum protocols. Quantum fidelity 
is the standard tool for assessing the success of quantum tele¬ 
portation IMl , where an unknown state is destroyed in one 
location and reconstructed in another (^e Ref. 1113] for a re¬ 
cent review). In quantum cloning where an unknown 

state is transformed into two or more (imperfect) clones, quan¬ 
tum fidelity is the basic tool to quantify the performance of a 
quantum cloning machine. Quantum fidelity plays a central 
role in quantum metrology IllSl fl^ . where the goal is to find 
the optimal strategy to estimate a classical parameter encoded 
in a quantum state. Similarly, it is important in quantum hy¬ 
pothesis testing where the aim is to optimize the dis¬ 

crimination of quantum hypotheses (states or channels). 

An important setting for all the above tasks is that of 
continuous-variable systems which are quantum sys¬ 

tems with infinite-dimensional Hilbert spaces, such as the 
bosonic modes of the electromagnetic field, described by po¬ 
sition and momentum quadrature operators. For these sys¬ 
tems, Gaussian states 1123 are the most typical quantum states 
in theoretical studies and experimental implementations, so 
quantifying their similarity is of paramount importance. The 
derivation of a simple formula for the quantum fidelity be¬ 
tween two arbitrary bosonic Gaussian states is a long-standing 
open problem with a number of partial solutions accumu¬ 
lated over the years. We currently know the solutions for one 
mode ll25l - |y^ and two modes II28|] . A simple formula for 
multimode Gaussian states is only known in specific cases, 
namely when one of the two states is pure 1291] or for two 
thermal states 1^ . 

Here we solve this long-standing problem by deriving a 
computable formula for the quantum fidelity between two ar¬ 
bitrary multimode Gaussian states which is simply expressed 
in terms of their first- and second-order statistical moments. 
A key step for this derivation relies on the adoption of an 


exponential Gibbs-like representation for the Gaussian states, 
which has been used recently to evaluate the fidelity between 
fermionic Gaussian states and which allows us to sim¬ 
plify many calculations. We also provide a recipe for express¬ 
ing the quantum fidelity in terms of symplectic invariants, 
showing specific examples with one, two and three modes. 
The new formula for the fidelity allows us to easily derive the 
Bures metric for Gaussian states, therefore generalizing quan¬ 
tum metrology to multimode Gaussian resources. Similarly, 
we discuss how quantum hypothesis testing can be extended 
beyond two-mode Gaussian states. 


II. PRELIMINARY NOTIONS 

Consider n bosonic modes described by quadrature oper¬ 
ators Q - {x\,... ,Xn,Pu. ■ ■ ,PnY, satisfying the canonical 
commutation relations 

[e,e^] = /Q, (1) 

where 1 is the n x n identity matrix. The coordinate trans¬ 
formations Q' - SQ which preserve the above commutation 
relations form the symplectic group, i.e. the group of real ma¬ 
trices such that SQ.S^ - Q. 133 . 

Let us denote by p an unnormalized density operator of 
the n bosonic modes. Its normalized version is denoted by 
p = p/Zp, with Zp = Trp being the normalization fac¬ 
tor. For a Gaussian state M23 . the density operator p has a 
one-to-one correspondence with the first- and second-order 
statistical moments of the state. These are the mean value 
M := {Q)p - Tr(2p) e and the covariance matrix (CM) V, 
with generic element 

V/il - j({Qk - U/c, Ql - M/))p , (2) 

where {,) is the anticommutator. Equivalently, we may use the 
following modified version of the CM 

W := -IViQ.. 


(3) 
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According to Williamson’s theorem, there exists a symplec- 
tic matrix S such that 

V^S(D®D)S^, D = diag(vi,...,v„), (4) 


where the symplectic eigenvalues satisfy Vk > 1/2. Corre¬ 
spondingly, the matrix W transforms as and its stan¬ 

dard eigenvalues are +Wi- where Wk - 2vk > 1 ■ 

In AppendixiAl we show that an arbitrary multimode Gaus¬ 
sian state with mean u and CM V can be written in the expo¬ 
nential form 


p - exp 


1 


(2 - uYG{Q - u) 


Zp = det IV 



where the Gibbs matrix G is related to the CM by the formulae 


G - 2 /Q coth *( 2 y/Q), V - - coth( -^ | /Q. 


Equivalently, we may consider the following relations 


jnc 


W-t 

wTI’ 




t + e 


iSlC 


1 - ’ 


( 6 ) 


(7) 


we use the notation AjB when A and B commute - 

see Appendix iBl for more details. Although the matrix G is 
singular for pure states (so one has to deal carefully with this 
limit), the introduction of the representation in Eq. (IS) signif¬ 
icantly simplifies the calculations, and all the final formulae 
are valid in general, i.e., for both mixed and pure states. 


III. FIDELITY FOR MULTIMODE GAUSSIAN STATES 


The quantum fidelity between two arbitrary quantum states, 
pi = pi jZp^ and p 2 = pi/Zp,, is given by 


Tipupi) := Tr^-y/ -v^p2 


( 8 ) 


where ptot := ■\fpip 2 2 jpi- We consider two Gaussian states, 
pi with CM Vi and mean ui, and p 2 with CM V 2 and mean 
U 2 . The Gibbs matrices Gi and G 2 are readily obtained from 
Eqs. (|6]l and 0. The advantage of the Gibbs representa¬ 
tion 0 for the calculation of the hdelity is twofold: hrstly, 
it makes the evaluation of the operator square root in Eq.0 
straightforward, and secondly, one can use the algebra of 
quadratic operators to hnd ptot in a closed form. 

As we show in Appendix O given two generally-displaced 
Gaussian states, the formula for their quantum hdelity can be 
directly expressed in terms of i5„ := U 2 - ui and their CMs, Vi 
and 1^2. In fact, we hnd 


T(PI,P2) = To(yi,y2)exp 


1 ^ 


V2r6u 


, (9) 


where the term 'FoiVi, V 2 ) depends only on Vi and y 2 and is 
easily computable from one of the two auxiliary matrices 


More precisely, we hnd 


Wl,V2) = 




Plot = det 


= det 


1 + (^auxO)-^ ^ 




Va, 


( 12 ) 

(13) 

(14) 


Note that the asymmetry of y^ux and Waux upon exchanging 
the two states is only apparent and comes from the apparent 
asymmetry in the dehnition of Eq. 0. One can check that 
the eigenvalues of Eaux and Waux. and thus the determinants in 
Eqs. (fTSl) and (fl4l) . are invariant under exchange. 

We remark that the formula of Eq. 0 is valid for arbitrary 
(generally-mixed) multimode Gaussian states with arbitrary 
hrst- and second-order moments. In the specihc case where 
one of the states is pure (say pi), we have Vi - 1/2 which 
implies Eaux = 1/2 and Etot = 1, therefore recovering the 
recent result of Ref. (in different notation 


IV. EIDELITY IN TERMS OE SYMPLECTIC INVARIANTS 


The hdelity can be expressed in terms of symplectic invari¬ 
ants associated with the second-order moments of the Gaus¬ 
sian states. Consider the notation with the W-matrices, so that 
Etot is given by Eq. (fT4l i. The standard eigenvalues of Waux are 
+w™’‘, where > 1 JUt]. As a consequence, we may write 


Etot 


fl [<“ + 


(15) 


Thus, the problem reduces to hnding the eigenvalues of VEaux- 
Eor this, let us consider the characteristic polynomial 


^(d) = det(Tl-iyaux), ( 16 ) 


which is clearly a symplectic invariant since VEaux transforms 
as ^VEaux^ * under symplectic transformations. Usin g th e 
identity dete^ = and the Cayley-Hamilton theorem 11^ . 
we may write ;^f(/^) as a polynomial function of 

fork = (17) 

which are also symplectic invariants with 4 > Ij for k > j. 
Thus, for n modes, we can compute the n invariants Gk and 
subsequently solve the polynomial equation ;^f(/^) = 0, whose 
roots are the eigenvalues to be used in Eq. ( fTSl ). 

Note that the invariants Gk can be connected with other in¬ 
variants. Eor instance, one can easily check that 

^(0) = (-1)"^, ;^(1) = (-1)«^, (18) 

where A := det(yi -H y 2 ), E := 2^" det(GyiQy 2 - 1 /4) and 


Vaux = n^(Vl + V2)^‘ (J + V 2 SIV^ , 

Waux ;= -2yaux/o = -{Wi + W2r\i + iy2 Wi). 


( 10 ) 

( 11 ) 


A := 2^" det(yi -H /Q/2) det(y 2 -H /Q/2) (19) 

are the invariants considered by Ref. ll^ . Using Eq. ( fTSl ), one 
can easily express I 2 and 4 in terms of E, A and A. 
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V. EXAMPLES 


Let us show some examples with n - I, 2 and 3 modes. 
For single-mode Gaussian states, we derive ;t'(/l) - - I 2 I 2 , 

so that = V.^ 2 / 2 . Equivalently, we may compute I 2 I 2 . - 
1 -H A/A so that we retrieve the known result 112514^ 




1 

VAh-A- Va 


( 20 ) 


For two-mode Gaussian states, we derive = (/j - 2/4 - 
4 / 2 /!^ -H 8/l^)/8 with solutions 

wr = \ ^/2 ± - 1\. ( 21 ) 

Once plugged into Eq. (fTsT l. we have the hdelity in terms of I 2 
and l\. The latter invariants can then be expressed in terms of 
r/A and A/A, so that we retrieve the known result 

, y 2 ) =- / ^ =. ( 22 ) 

Vf-HVA-y(Vf-HVA) - A 


moments u and V, and p 2 - p + dp, with statistical moments 
u + du and V + dV. Then, the Bures metric is given by 


ds^ = 2[1 -T(p,pH-iip)] 


du^V ^du 
4 



(26) 


where 6 ATx[dV{A£.Y + -Ca) ^dV], X.aX AXA, and the 
inverse of the superoperator AX.v + -Cn refers to the pseudo¬ 
inverse (see Appendix|E|for the proof). Note that a result 
equivalent to Eq. (|2^ has been derived in Ref. using a 
different method based on the computation of the symmetric 
logarithmic derivative. 

Numerically, the easiest way of evaluating the inverse of the 
superoperator in 6 is using the VT-matrices and performing the 
calculations in the basis in which W is diagonal. In the basis 
where W is diagonal, then 




u 


dWijdWji 
WiWj - 1 


(27) 


and the sum is taken over the elements such that w,wj 4^ 1 . 
For pure states, we simply have dpure = Tr(y“' dV y^' dV). 


B. Multimode quantum metrology 


For three-mode Gaussian states, the characteristic polyno¬ 
mial may be written x - + pt + q, where 


f = 42 - I2I6, p = ^ 


h hl4_h 

4 ’ ^ 108 12 6 ■ 


(23) 


The solutions of the characteristic equation = 0 are real (see 
Appendix iDli and given by 


72^, / - 
= -i- 1 - 2 .,/— cos 


e-2n(k- 1) 


(24) 


where 0 ;= arccos |^3 yf3q{2p yj-p) and k — 1 , 2 ,3 (in par¬ 
ticular, note that w™’' = V 72/6 for p - 0). To the best of 
our knowledge, Eqs. (|2^ and (l24l) . together with Eqs. (01 
and (flSl l, provide the hrst expression for the quantum hdelity 
between two arbitrary three-mode Gaussian states. 


VI. IMMEDIATE IMPLICATIONS 
A. Geometry of Gaussian states 

Once the quantum hdelity is expressed in terms of the hrst 
two statistical moments, we can easily compute the Bures dis¬ 
tance between two arbitrary multimode Gaussian states, pi 
and p 2 , which is given by 


Let us consider a real parameter 6 which is encoded in a 
multimode Gaussian state pg. To estimate 6 with high pre¬ 
cision, it is necessary to distinguish the two inhnitesimally- 
close states pg and pg+de for an inhnitesimal change dO. As¬ 
sume that N copies of the state pg are available to an observer, 
who performs N independent measurements to obtain an unbi¬ 
ased estimator 6 for parameter 9. Then, the mean-square error 
affecting the parameter estimation Var(0) {{0-0)^) satishes 
the quantum Cramer-Rao (QCR) bound Var( 6 () > [A//(0)] ', 
where H(6) is the quantum Fisher information (QFI) ifT^ . The 
latter can be computed from the hdelity as 


H{9)^ 


8 [1 - ‘Fipg, pg+de)] 

dfi 


(28) 


Thus, for any parametrization of the Gaussian states, we can 
easily compute the hdelity 'Fipg, pg+de) using Eq. (0 and, 
therefore, the QFI in Eq. (l28l l. 

More generally, suppose that the Gaussian state is labelled 
by a vectorial parameter with m real components, i.e., 6 = (ft) 
for / = 1,..., m. In this case, the performance of the parame¬ 
ter estimation is expressed by the classical covariance matrix 
Co\ij{9) ;= (OiOj) - (9i)(9i), which satishes the matrix ver¬ 
sion of the QCR bound Cllll Cov( 6 l) > [A//(0)]-'. Here 
the QFI is a matrix with elements HpiO), which can be eval¬ 
uated from the Bures metric. In fact, for any parametrization, 
we may write Eq. (|26| as ds^ - gij(9)d6id6j and show that 
Hij{0) = Agij{0). 


£>B(pi,P 2 ) = 2[l-r(pi,p2)]. (25) 


C. Multimode quantum hypothesis testing 


Form this expression we can derive the Bures metric 
by expanding the hdelity. In fact, let us consider two 
inhnitesimally-close Gaussian states pi = p, with statistical 


An efhcient computation of the quantum hdelity is cru¬ 
cial for solving problems of binary quantum hypothesis test¬ 
ing ll^l2ll] with multimode Gaussian states. These problems 
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may occur in the basic scenario of quantum state discrimi¬ 
nation, where two Gaussian states must be optimally distin¬ 
guished, or in the setting of quantum channel discrimination, 
where two Gaussian channels must be distinguished by as¬ 
suming Gaussian sources and input energy constraints. In par¬ 
ticular, the latter formulation is very important in a variety of 
quantum technology protocols, such as remote quantum sens¬ 
ing of targets, i.e., quantum illumination ImlM^ . and quan - 
tum reading of classical data from optical memories 

Consider N copies of two multimode Gaussian states, 
and p®^, with the same a priori probability. The minimum 
error probability Perr{N) in their statistical discrimination is 
provided by the Helstrom bound ifs^ . which is typically hard 
to compute for mixed states. For this reason, one resorts 
to other computable bounds, such as the quantum Chernoff 
bound |0-S2l or fidelity-based bounds Thanks to 

our result the latter are now the simplest to compute. 

For any number of copies N, we may write 

1 - Vi - .r«s,.A)i» 

-2- - Ferr(A^) < -^ 

In particular, the lower bound in Eq. ( |29] ) is the tightest known. 
Note that Eq. (l29l) can be derived by using the known result 
for single copy (N - 1) and then applying the multiplica¬ 
tive property of the fidelity under tensor products of density 
operators, so that T(pf^,P 2 ^) - T{pi,p 2 )^■ 

The computation of the quantum fidelity is also important 
for asymmetric quantum hypothesis testing where the two 
quantum hypotheses have unbalanced Bayesian costs In 


this context, the quantum fidelity can be used to estimate the 
quantum Hoeffding bound which quantifies the optimal 
error-exponent associated with the rate of false negatives. 

VII. CONCLUSIONS 

In this work we have solved a long-standing open prob¬ 
lem in continuous variable quantum information by deriving a 
simple computable formula for the quantum fidelity between 
two arbitrary multimode Gaussian states. Our main formula 
is expressed in terms of the statistical moments of the Gaus¬ 
sian states, but another formulation is also given in terms of 
suitable symplectic invariants. By using our formula, one 
can extend the study of quantum teleportation, cloning, quan¬ 
tum metrology and hypothesis testing well beyond the stan¬ 
dard case of two-mode Gaussian states to consider multimode 
Gaussian resources, with unexplored implications for all these 
basic quantum information protocols. 
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Appendices 

Appendix A: Exponential formula for Gaussian states 

Here we show the formulae in Eqs. ©-(Ell. The hrst step 
is to introduce the symplectic action of a real function / on a 
CM and how it can be computed in terms of standard matrix 
functions when / is odd. After this preliminary step, we start 


by noting that, for thermal states (having V - D® D), we can 
easily write Eq. © with u - 0 and 

G^g{D)®g{D), g(v) = 2coth-‘(2v). (Al) 

Then, we generalize the formula to zero-mean Gaussian states 
with arbitrary CMs by noting that QGQ transforms as V under 
symplectic coordinate transformations Q’ - SQ. This prop¬ 
erty allows us to use the symplectic action gt(v) which leads to 
Eq. ©. Finally, we include displacements to extend the result 
to arbitrary mean values and we compute the normalization 
factor. 


1. Symplectic action and its computation 

Then, let / : R ^ R be a function. The symplectic action 
/* on the CM V is dehned by il^ 

MV)^S[f(D)®f(D)]S\ (A2) 

where /(D) = diag[/(vi),/(v 2 ),... ,/(v„)] acts as a standard 
matrix function. Here we prove that, if / is an odd function 
f(-x) = -fix), then 

MV) = /(y/G)/Q. (A3) 

Let us start by proving that Eq. (IA3I) satisfies the identity 
MSVS^)^SMV)S^ . (A4) 

In fact, we have 

MSVS^) = fiSVS'^in)iQ. = fiSViQ.SMiVl 
= SfiViD.)S-^in = 5/(y/Q)/G5^ 
^SMV)S^, 

where we use the basic property f{S y5“') = SfiV)S~^. 

Because of Eq. (IA4b . without loss of generality, we can fo¬ 
cus on the case where V is in diagonal Williamson form, i.e., 

V - D®D, D - diag(vi, V 2 ,..., v„), 

and we assume that v, vj for i + j. One can easily check 
that the matrix 

y = (D © D)iQ. 

is Hermitian, so it can be cast into the diagonal form by a 
unitary matrix U. It turns out that U is independent on v, and 

V^U\D®-D)U, (A5) 

with eigenvalues ±v,. If / is an odd function, then 

/(y) = M[fiD)®fi-D)]U = M[f{D)®-fiD)]U. 

The latter matrix has the same structure of V in Eq. ( IA5I ). Be¬ 
cause U is independent on the diagonal elements, then 

/(y) = [/(£>) e/(D)])Q, 

which gives 

/[(D©D)/Q]/Q = /(D)©/(D). 

This is Eq. (IA31 l up to a symplectic transformation S. 
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2. Proof of the exponential formula 


Let us now show that the Gibbs exponential formula of 
Eq. 0 can describe an arbitrary Gaussian state (not just a 
thermal state). We start by considering a single-mode thermal 
state p - In this case, we can write 






1 dz 


1 


Zdg e« - 1 


(A6) 


i.e., matrices V and QGQ transform in the same way un¬ 
der symplectic coordinate transformations. As a result, they 
can be related by the symplectic action of the function in 
Ea. (lAT2l l. 

In fact, for thermal states, we may write 

y = D©D, (QGQ) =-g(D©D). 

Then, for an arbitrary symplectic transformation S , we have 


In our notation, a - (x + ip)I V2 so that a^a = — 5 and 

v(g):=<x2) = </) = (a1a) + i/2. (A7) 

Therefore, from Eqs. (IA 6 I 1 and (IA71 i. we derive 


Kg) = 5Coth|. 


(A8) 


Thermal Arbitrary 

D®D 4 y = 5(D©D)5^ 

-g(D®D) 4 QGQ = 5 [-g(D©D)]5^ =-g*(V). 

Thus, using the symplectic action g,, dehned from Eq. (IA12I) . 
and its inverse v», dehned from Eq. (IA8I) . we can derive the 
relations 


In terms of the quadratures, the thermal state reads 

p = (A9) 

and its normalization is given by 

Zp = Zpe-i - , ^ z(g). (AlO) 

e2 — e ^ 

Note that the purity is given by 

Trp^ = Zp 2 /Zj = z(2g)z"^(g) = tanh(g/2) = ^Kg)”', 

so that the vacuum corresponds to g —> cx) or v —» 1/2. 

The previous representation of Eq. ( IA9I ) can be generalized 
to a multimode thermal state of n > 1 bosonic modes. This 
state has its CM already in the diagonal Williamson form 

V-D®D, D = diag(vi,..., v„). 

Thanks to the tensor product structure, we can write 


Here G := diag(gi,... ,g„;gi,... ,g„), where the diagonal el¬ 
ements are given by g, = g(v,), where 

g(v) = 2coth^'(2y) (A12) 

is the inverse of the function in Eq. (IA8b . Compactly, we set 

G = g(D)©g(D). 

Now, we study how G and V transform under coordinate 
transformations Q' - SQ. We have V - SVS^ and 

G' = = QgQGQ5^Q, (A13) 

where Eq. (IA13I I comes from imposing Q^GQ - Q'^G'Q' in 
Eq. (lAl lb . Erom Eq. (IA13b . we see that 


G = -2Q coth.,^(2y)Q = 2/Q coth '(2y/Q), 


and 



QGQ\ 



/Q, 


where we also exploit Eq. (IA3b . These formulae correspond to 
those in Eq. (|6]l given in the main text. The additional formula 
in Eq. O is obtained by considering that W - —2ViD.. 


a. Extension to non-zero mean 

The next step is to include the presence of a generally non¬ 
zero mean value in the exponential expression of Eq. (lAl lb . 
Eor an arbitrary u e R^", consider the displacement operator 

D{u) = 

which satishes D{u)^ - D{-u) and D{u)QD{u)^ - Q -i- u. 
By applying this operator to Eq. (lAl lb . we can generate an 
arbitrary Gaussian state with non-zero mean 

p = D(-u)e-'^Q"^<^D{u) = e-i(e-«Zc(G-«)^ 

This is easy to double check. Let us set 

p = D(-u)pcD(u), Pc 

Eirst note that Zp = Zp^. Then, we can verify that 

g-k(Q-ufC(Q-u) g-/G^CG 

Tr[2---] = TiiQDi-u) --- D{u)] = 

e4e"CG e-hO^CQ 

Tr[D(M) 0 D(-m)— -] = Tr[(e + u) — -] = «, 

Zp Zpj. 


(QGQ)^ 5(QGQ)5K 


i-e. {Q)p = M. Similarly, Vij = ^{{Qi - m/, Qj - Uj}}p. 
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b. Normalization factor 


The trace of an unnormalized Gaussian state p is written in 
Eq. (IA6b via the function z(g) - defined in 

Eq. (lAlOb . When G is diagonal (i.e. V is diagonal) then 

^P^W^^Sj)- (A14) 

j 

Now we write Eq. (IA14b in a coordinate independent form. A 
generic G can be obtained from a diagonal G via a symplectic 
coordinate transformation, because of the property (IA2b of the 
symplectic action, and because det 5=1, one has 

Zp = Vdet[z.(G)] = det[z(G /Q) /G] ^ det[z(/Q G) /Q] 

= det[(e'“'^''2 - . (A15) 

Moreover, z(g(v)) = - j. It is simple to prove that 

Zp=n =n i 

j j 

= det[ydiag + /G/2]*/2 (A17) 


where Vdiag = diag(vi,..., v„; vi,..., v„). Since a general V 
can be written as E = 5 Ejiag S ^ and det 5 = 1, then 


Zo = detly + 


1/2 


where we used the fact that 5 Q.S ^ = G. By replacing W - 
-2y/Q, we also get 


Zp = det 


11 

2iQ j 


Appendix B: Computations with Gaussian states 
1. Product of two Gaussian states with zero mean 


Although the product of two Gaussian states can be read¬ 
ily evaluated thanks to the result of lHH, in this section we 
provide a self-consistent proof. 

By using the Baker-Campbell-Hausdorff identity, we can 
write the product of two zero-mean Gaussian states as 


The above identity is a consequence of the algebra 


Q^GQ Q^G'Q 
2 ’ 2 


^e^(GQG' 

Q^G"Q 
2 ’ 


g'qg) q 

(B2) 


where 


G = -i'G/, & = -iQ/', G" = -inJ" 


and J" = [7, J']. Because of the above identity, we can write 
the Eq. (IB lb with e'' - , namely 




(B3) 


Now we can express the composition rule of Eq. (IB3b in 
terms of the CMs V and V' of the two states p and p'. Erom 
Eq. (©, we have 

1 H- 1 ^ ,or -2yiQ - 1 

2 e'£2G _ 1 -2y/G -H 1 

In terms of IT = -2y/Q, W and W" , we may write 


ly" 


ginc'g/nc + I 


iv'-ll ^ w+1 
w'+H w-11 


ginG'g/nc _ I w'-l _ w+l 
w'+i w-i 

1 - 2(iy' H- 1)-' + 1+ 2(W - 1)-^ 
~ 1 - 2(iy' H- 1)-' - 1 - 2(iy - 1)-' 

1 -(ly' H-1)-' -H(iy- 1)-' 

(ly' H-1)-* -H(iy-1)-' 

^ ^ t-2{w' + ir^ 

^ (W + 1)-' +{W- !)-'■ 


In the above equations, we fix the notation 4 = AB ' when 
[A, B] + 0. Using the Woodbury identity 

(A -H = A”' - A^'(A^^ -H 


we derive 

ly" = 1 -H (w' + 1 - 2 ii)[i - (w' + iy)-‘(iy' H- 1 )]. 


Then, using another straightforward matrix equation 
(A -H B)^‘A = 1 - (A -h B)"'B, 
with A = ly' -H 1 and B — W - 1, we find 

ly" = 1 iw' - t)(W' + wr\w - 1). 


Therefore 

y" = 


iQ. 

y 



(y' -H y)-' 



(B4) 

(B5) 


(B6) 


Note that the squared of a Gaussian state p^ has - 
2G and its CM can be computed directly from the previous 
Eqs. (IB5b and (IB6b by setting W — W and y = y'. It is easy 
to check that we get 

^(2)^ 1 iy(2) ^ 1 


2. Square root of Gaussian states 


Given a Gaussian state p, its square-root yfp is a state with 
G ^ G/2. The CM Esq of ^Jp can be written in terms of the 
CM y of p by concatenating functions 


Vsq(v) := v(g(v)/2) = 



V. 


(B7) 
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Notice that, because v > 1 /2 one might be tempted to sim¬ 
plify Vsq(v) into the expression x + - 1 /2. However, the 

latter function is not odd, so it produces wrong results when it 
is used for the symplectic action. Eq. (IB71 i is the correct one. 
When V is Williamson-diagonal, so it is Vsq and the diagonal 
elements are given by [Vsq];; = Vsq(v;). Since V and Vsq trans¬ 
form in the same way under symplectic transformations, for 
any general (non-diagonal) V, the relation between V and Vsq 
can be obtained with the symplectic action 


Esq — Vsq^*(E) — 



(yn)-2 

4 


-H 1 


V . 


By replacing W - -2y/Q, we finally derive 

Wsq = (Vl - W-2 -H l) W. (B8) 


3. Extending the product formula to Gaussian states with 
non-zero mean 


Then, using the above result 

j(Q-uV 
- p(G, u) 




(B12) 


where 


1 


- ^pTt„Gin _ _ JQC\t^-iQ.G _ |^-1 ^ 




= -r iQ.{e 


]l)-‘/Q(e“ 


T;n,tJO.G _ I, \-\i„-iSlG _ „inC\[„-iQ.G _ |^-1 ^ 




= - 1)^' - - iy^){ 


= -r iQ. 


W+l W-l 


1 r 


When an operator linear in terms of Q is introduced, the 
algebra in Eq. (IB2t has to be extended. It turns out that 




i£lGQ , 


[uQ, Q\ - wiQ.v . 


(B9) 


Therefore, D(m)2D(m)^ = Q+u, and using Eqs. (IB2I) and (IB9I) . 
we may write the identities 

D{u)D(v) = D{u + 

e^Q^Gag^iQ^GQ ^ ^inGg (gio) 


4. Decomposition of displaced Gaussian states 

Using the previous identities we may write 

p = e-i(G-HEG(e-«) ^ g-H^'/nCg-jG’'ceg/H’'/nG 
^ g-«''inGg«''/nc'“°Gg-|G''cG 
^ g-H’'inGg«’'e“‘'inGg-iG’'cQ 
^ ^-u'^inQ^{e-“‘°ufinQ^-\Q'^GQ 
_ g 1ine-'“°Mg(e-'“°M-M)’'/nGg-1 G^'gg 

Let { - - u, i.e. 

u = 

Note that 

u - ‘%[iyu] 

= - e‘^^)u . (B11) 


Appendix C: Proof of Eq. @ 


We start by considering the undisplaced case where ui - 
U 2 - 0. This assumption will be relaxed in Annendix ICSI 
The total state ptot := -\/pi has CM Vtot (Wtot) and its 

Gibbs matrix Gtot can be derived by applying the composition 
rule of Eq. (IB31 l and noting that sjp has G/2. Thus, we have 


exp (/QGtot) = exp 



exp (/QG 2 ) exp 



(Cl) 


Using the expression of the partition function Zp in Eq. (|5]), 
the relation between the CM V and the Gibbs matrix in Eq. (|6l) 
into TipuPi) = Z^J yjZp^Zp^, we may write 


TipirpiT^ 


dct _^“^^^ 101 / 4 ^ dct _^“^^Giot/4^ 

^^/nG2/2 — 

det(e'“‘^“^2 _ ]ljdet(e'“‘^“/2 _ 
det(e'^Gi _ t)det(e'^G2 _ |) 

det (jgf ^g/nc .„,/2 _ 

det(e'^Gi _ t)det(e'^G2 _ |) det(e'^G,ot/2 + j) 
det^e'^^2 — det^e'^‘^“‘^2 _ 

det(l - e“'^Gi)det(e'^G 2 _ 1 ) det(e'^Gioi /2 + j) 
[r(Gi,G2)Etot]-4 (C2) 


where 


r(Gi,G2):= 


|det(]l -e '^Gi)det(e '^‘^2 _ |) 

det (e'^ G 2 _ g-ffi Gi ^ det/Q 


Efnt 


/det(e'^G,oi/2 + 
det (e'^Giot/2 — 1) 


det iQ.. 


(C3) 
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Now it is easy to check that 

r(Gi,G2) = 


-^det[yi + V’ 2 ] 


(C4) 


By contrast, the computation of Ftot is more difficult. Using 
Eq. O we may write ^’tot in terms of fktot as follows 


Ftot = det 


(V 


1 - Wj + 1 


or, equivalently, in terms of Vtot as follows 


= det 


1 + 




-2 


Vtr 


1/4 


1/4 


(C5) 


(C6) 


Let us compute Vktot as a function of Wi and For this 
we iterate the composition rule in Eq. (IB51 l and we use the 
following relations for the VE-matrix of the square-root state 


(C7) 

(C8) 


lUsq = (V]1-VE-2 -h1)vE, 

Let us start by applying Eq. (IB5I) twice. We have 
W" ^1+ (W 2 - ll)(Wlsq -H W2r‘(Wlsq - 1), 

Wtot = 1 + (Wlsq - l)(Wl,q + W"r\W" - 1) 

= 1 IXWlsq-H W"r‘x 

(W 2 - i)(Wuo + W2r\Wu. - 1). 


Now the next step is to apply the Woodbury identity and (A * - 
= A(A + B) '‘B multiple times, so that we have 

[Wlsq + 1 + (W 2 - + W2r'(Wl,q - 1)]-' = 




(Wlsq - 1) 

and we may write 

Wtot = 1 + 


(Wl,q + 1^2)-' + 


1 Wlsq + 1 


W 2 - 1 Wlsq - 1 


-1 


Because Wisq-HW 2 = Wisq-HW 2 -HWi-Wi and IdUisq-W-^q) = 
VEisq - VEi, we may write 


Wtot = Wl - (Wlsq - Wl)(Wi + W2r‘(W^lsq - Wl). 

This is already a simple expression, but it can be further 
simplified. Let us write its inverse 


= 

''tot 


1 


Wi 


1 


Wlsq - Wl \ (Wlsq - Wl )2 Wl + W 2 / Wlsq - Wl ' 

Using Eq. (IC7l l we may write 

(Wlsq- Wl)2 


Wi 


= Wi - Wf', 


which, replaced in the previous expression of leads to 
W^lsq - Wl ,-l 


w~ = 


X (Wi + W 2 ) 


(VEi - Wf' - VEi - IE 2 ) 
1 


Wlsq - Wl 


= -(Wisq - Wi)(i -H W2WirHwi + W 2 ) 


Wlsq - Wl 


= (Wl.q - Wl)W,- 


1 


Wlsq - Wl ’ 


where 


W = -- 

aux 


Wl 4-VE 2 


(1 -HIE 2 IE 1 ). 


(C9) 


(CIO) 


Because in Eq. (IC51 l there is a determinant of matrix func¬ 
tion, such expression is invariant under MVE/W ' transforma¬ 
tions (with non-singular M). Therefore, we can use Waux in 
the place of VEtot in Eq. (IC5b . In other words, we may write 


(W 2 -II)-', 


Etot = det 






1/4 


I_ 

+ 

= det 

2 



(Wlsq -H W 2 ) -H- 

*'*‘1 1^2-1 Wlsq - 1 . 

X 


[S 4 J 



1/4 


(Cll) 

(C12) 


(Wlsq + W2)-XH^lsq - 1) 

= fVlsq - (fVlsq + fV2)2f^‘(fVlsq - 1) , 


where we have used Waux = -2yauxif^- Combining Eqs. (IC2I) . 
(IC4I) and (IC121 i. we obtain Eq. (|9|l, (fOT l and (fT4l i. 


where 


X = Wlsq -H iy2 -H 


Wlsq - 1 
Wlsq + 1 


(iy2 -1) 


1 


Wlsq + 1 
Wlsq 
Wlsq + 1 


(1 + + 2iyisqiy2) 

, 2iyisq 

-1 , w. , ow.x _ q + ^ 2 ), 


(^Isq + ^Isq + 2iy2) = 


Wlsq + 1 


and we have used Eq. (IC8I) . Therefore 


Wto, = Wisq - -(Wisq + iy2)(iyi + iy2)-'(iyisq - wr/q). 


1 . Comment for pure states 

The most important result of the previous sections is the 
similarity transformation which relates VEtot and Waux: 

Wtot = (Wisq - (C13) 

VEl.sq - Wl 

However, whenpi is pure VEisq = VEi so the above transforma¬ 
tion is singular. The purpose of this section is to show that the 
final result ( IC12l i is consistent even when the matrix VEisq - VEi 
is singular. 
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To simplify the notation we assume that pi is a pure state, so 
the symplectic eigenvalues v' are equal vj = 1 /2, Vi, although 
the following argument can be easily generalized to the case in 
which only few eigenvalues are equal to 1 /2. Because Eg. (ICS b 
is basis independent, we perform the calculation in the basis 
where Wi is diagonal and we write 


Wi(e) ^ eDu Di=lle(-ll). (C14) 

6-^1 

Since Ea. dCSb depends only on the eigenvalues of Wtot and the 
eigenvalues are smooth under perturbations we can write 


Etot = lim det 




1/4 


(C15) 


where lTtot(6) refers to VTtot with Wi substituted by Wi(e). Eor 
any e < 1, it is VTisq - Wi - Vl - e~^Di so the similarity 
transform (1C 13b is well defined and (1C 15b can be replaced by 
(ICl lb . Although the matrix VTisq - Wi is singular for e ^ 1 
its dependence cancels out, while VTaux is well-defined even in 
the limit e ^ 1. 

This is confirmed by the fact that (IC5b reproduces the 
known results when pi is pure. In the next section we 
expand this point to simplify the numerical treatment of the 
singular case. 


2. Treatment of the singular case 


In this section we devise a strategy that helps the numerical 
treatment of the singular case, i.e. when one or more sym¬ 
plectic eigenvalues of Vi and/or V 2 are equal to 1/2. Because 
the eigenvalues of ITaux are invariant under the exchange of 
the states pi p 2 , without loss of generality we assume that 
Vl is the state with the highest number of eigenvalues equal to 
1 /2. Let r be the number of pairs of symplectic eigenvalues 
of Vl equal to 1 /2. Since Vaux transforms under symplectic 
transformations, without loss of generality we can perform the 
calculations in the coordinate system where Vi is diagonal. 

Moreover, to simplify the notation, in this section we re¬ 
shape the matrices so that Q = ©y ^j. Therefore, we can 
can write Vi and V 2 in the block form where 


El 


_(i2r/2 0 

\ 0 D 


V 2 = 


A C 
B 


Q. = 


CO 0 

0 Cb 


(C16) 


where l 2 r, A, to are 2rx2r matrices, C is a 2rx2{n-r) matrix, 
and D, B, at are 2(n - r) x 2{n - r) matrices, D is diagonal 
with diagonal entries greater then 1 /2. Thanks to this block 
structure, with a long but straightforward calculation we find 


Va, 


/ 112./2 (5\ 

I 0 5/’ 


(C17) 


where the matrices C and B depend on A, B, C, D. Because 
of the block structure of Eq. ( IC17b . it is clear that VTaux has 
r eigenvalues equal to 1 and r eigenvalues equal to -1. In 
view of Eq. (ffSl) . these eigenvalues do not contribute to Etot 
and can thus be discarded. On the other hand, the eigenvalues 
^aux ^ | 3 g found by diagonalizing ITaux = -2iBcb. 

With a similar argument, hk -2r + 


3. Alternative Formula 


Note that in the proof of Sec. Owe can exploit the fact that 
det[/(y)] = det[/(t/yt/“')] for some invertible matrix U. By 
using Eq. (IClb into Eq. ( IC3b . we get 


det 1^ Ve'ti Cl /Zgin G 2 gin Ci /2 q. j j 

“ det[Ve'“'=>/V“G 2 e/nGi/ 2 _ 

and with either U = e'^Gi /2 or f/ = 

det 1^ j j 


det iQ, 


- 

^ tot 


det [ j j 

det 1^ ginG 2 q. £ j 
det y - 1 j 

Einally, after simple algebra, we may write 


det /Q 


det /Q. 


E,ot = det 


= det 






Ti 2 

1^21 


being 


/Q 


iQ. 


yi2 - —— +1 El -H — I (El + y2) * I y2 + 


/Q 


(CIS) 


(C19) 

(C20) 


and y 2 i = y] 2 . Note that, contrary to matrix Eaux, the new 
matrix yi 2 is not real. Because of the above derivation, 

Wtot = ^ g-^nG,/2y^22^e/nGl/2 

the matrices Wtot, Wn and Vy^i are similar. 

The relation between yi 2 and Taux is easy to obtain using 
the W matrices and applying the Woodbury identity. We find 


Wi2^t + {Wi - ]i)(iyi H- W2)'‘(iy2 - 1) 

1 

= (Wi - 1 ) 


1 


W 2 W 1 -Wi-W2 + i Wi H- iy2 

= (iy 2 - i)-'(iy 2 iyi + i)—l—(W2 - i) 

VTi + iy 2 


= (iVi - 1 ) 


Wi -Hiy 2 


(iy2iyi -H i){Wi - I)"' 


(iy 2 - 1 ) 


(C21) 


so that iyi 2 = -UW^uxU * for some invertible U, as we can 
see by comparing Eq. (IC21b with Eq. (IClOb . 


4. Exchanging pi and p 2 

The final result for the fidelity, Eq. dCl lb . depends on the 
matrix Waux which is not symmetric upon exchanging pi and 
p 2 . This is due to the apparent asymmetry in the definition 
of the fidelity ([8]). However, we show here that ( ICl lb is in¬ 
variant under such exchange, even though VTaux is not. In¬ 
deed, thanks to the results of the previous section, if F{W) = 
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det [(Vl - lV-2 + l) wf\ then F,,, = F{W,oi) = ^’(Waux) = 
■^( 1 ^ 12 ) = F{W 2 \)- Because IVaux is similar to Wn (apart from 
a global sign), which again is similar W 21 , if we exchange p\ 
and p 2 , the resulting VKaux (with indices 1 and 2 swapped) is 
similar to the original one. Therefore, (ICl II) is invariant under 
such exchange. 


By replacing the latter expression into Eq. (IC221 i. we derive 

^ g <5,,gI5„)’'/nWi2(5, 6.) 

- )«„g j<5,( in(l-e‘^°i )lVi2( )S„ 

The term /ftot can be simplified noting that 


5. Derivation of the fidelity for displaced Gaussian states 

Consider displaced Gaussian states, pi having Gibbs matrix 
Gi and mean value iq, and p 2 , having G 2 and iq. Then 


-inc, _ /nc. ^ m + 1 _ Wi-a . 

ITi - 1 ITi + 1 ■ 

(l-^>'“^')lTi2(l-e-'“^‘) = -7;T^ 

Wi + 1 


4Wi 

Wf-l' 


IE 12 


ITi - 1’ 


T(pl,P2) = 










V^^lot 


- F^iPCi^PGi)' 


VPlot 

7 _’ 

^ V^^lot 


where Tipc^ ^PGi) is the fidelity (already computed) between 
two undisplaced Gaussian states, i.e., with Gibbs matrices Gy 
and G 2 but zero mean values. Therefore, we only need to 
compute Z^JZ^^. If we write 


andlTi 2 = 1+ IT 1 -1-(ITi + ]l)(lTi+ W 2 )'‘(W^i -1), which is 
a consequence of the identity (IB4b . Therefore, we may write 

Ki,i = 6liQ.{Wi + W2r^6u = -^liVi + V2r'^u, 
and finally 


T(p\,P2) 


(det[yi + y 2 ])‘^" 


^ — g 2^tl Wtoi) Giot ((? Wlol)+^lot 


then 


Appendix D: Proof that the solutions for the three-mode case 
are real 


v !— - Z I _ 

Moreover, from the definition one can see that 

gtr.„. ^ ^Ao. ^ ^PiP2 _ 
^PGio, ^PGiPG2 


As written in ( |2^ . in the three-mode case the characteristic 
polynomial (fThl) can be written as;t' - p + pt+q. The equation 
X - Q has real solutions if p < 0 and fA + jll < 0, which 
is simple to prove. Indeed, calling +w™’' the eigenvalues 
(C22) of one finds that hn = 2[(wf’')2'' -r + (w”’')^"]. 

Hence 


For the numerator we may write 

Zp,p, = Tr[pip2] = TT[D{-Ul)pc^D(ul)D{-U2)pG2D{u2)] 

= Tr[D(M2 - ui)pg,D{ui - U2)pc2\ 

where the phase in Eq. (IB 10b vanishes after the twofold use. 
Then calling 

^2 ^1 ? 

and calling G 12 the matrix such that 

g-inci2 _ g-incig-inc2 


one has 


Zp.p,=Tr[/^“V,^-^"'“V2] 

= Tr[/"'“eg-(.— 

_ S,yinQ^\sline-'^ 

Now, by using Eq. (IB12b we find 

^PlPl ~ ^ 


■15, 


'PC,2]- 


'''"’Tr[pG.PG2]- 


p = -^[(wr)"+(wr)"+(wr)" 

- (wf’‘)2(w^“)2 - (w™’‘)2(w^“)2 - (w^“)2(w^“)2] 

- («“)" - (wr)")'+(«“)' - (wr)T 

+ ((W^)^ - (W™)^)^ ] < 0. (Dl) 

Similarly, 

((w^“)2-(w^“)2)^] < 0 . (D2) 

Hence, the eigenvalues of IT™’' are real. The real solutions of 
= 0 are given by (l24l) . 


Appendix E: Derivation of the Bures metric 

Let us consider two infinitesimally-close Gaussian states 
Pi — p and p 2 - p + dp. The first state is parametrized by G 
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(or V) and u, while the second state is parametrized by G + dG 
(or V + dV) and u + du. Hence, up to the second order 


1 


1 1 


Vi + Vz 2y 1 + § 
and 


1 1 1 1 1 1 

- dV — + — dV — dV —, 

2y 2y 2y 2y 2y 2y 


dliVi + V2T^6 u ^ dJv^^dul2 . 
In a similar way, we find 
1 




t+W^+dWW 

2iy 1 


( 2 iy + dw) 


1 + iy2 
1 


1 + iy2 


dW 


W - 1 


ly 

-dW- - 


iy2 + 1 

iy2_ 1 


1 + iy2 1 + 1^2 I + ' 

Since the fidelity is an invariant, one can perform the calcula¬ 
tions in the basis in which ly is diagonal. Let us call ly the 
(diagonal) matrix W in this basis and dW the corresponding 
infinitesimal variation (non-diagonal). Then 


(^aux),.^. ^ 


2wi 


1 


TdWu- 


wl- 1 


' ‘ J 


Z 


1 


1 -I- w? 


dWikdWk 


Wk 


w]- 1 


1 -I- w2 1 H- w2 


To expand the expression 


■Finf - 


det(yi -H yz) 
det ITaux 


det(iy -H dWI2) 
one has to expand 


■ det 




W~^ + ]1 

aux “ ^ 




1 - + K' 


( 1 ), 


r(2) 


(El) 


(E2) 


in terms of the 0* order, first order and second order operators 
K^''\ Taking the square of Eq. (IE2t and calling 

^aux = 

the 2"“* order expansion of ly^ux^ we find the relations 

■^(0) = 1 “ yfo) 

K(\)K(0) + K(^o)K(i) = -y(i)y(0) - ^(0)12(1) 

K(2)K^o) + K(0)K^2) = -y(2)y(0) - V(0)y(Z) - yfi) - ■ 

These explicit calculation of is long but straightfor¬ 
ward. Once the operators K are known, from the expansion 

det(l + X) = gTrlogd+X) ^ ^Tr[X]-Tr[X^]/2 

of the three terms in Eq. (lEll i. we find 


Einf = exp 


1 vn dWijdWji 


i X 1 U VY i 

4 Zj TZ 


1 - WiWj 


i.e. 


T{p,p + dp) - exp 


8 16 Z-J 1 - vv, vv 

l] ' ■! ! 


The Bures metric is then given by 

1 dWiidWii 


ds^ - 


—du^ V ^ du -\— 
4 8 


8 WiWj - 1 

t / •' 


The above expression can be cast into a basis-independent 
form by defining the super-operator 


£aX = AXA . 


Indeed 


z 


dWijdWji 
WiWj - 1 


= Tr 


dW- 


XLw ~ 1 


-dW 


= Tr 


dW 


= -4Tr 


■Cw - 1 


dVQ. 


dW 


-Cw - 1 


(dVQ.) 


Using 


ZwidVQ) = -AVildVaVQ. = -A{J:v-LndV)Q ., 


we find 


Z 


dWijdWji 
4 WiWj - 1 


= 4Tr 

= 4Tr 

= 4Tr 

= 4Tr 

^2 _ 


dVQ. 


1 


dVZn 
dV 
dV 


A£.v-Ln + 1 
1 


A£,y£.ci + 1 
1 1 


■Cn A£vCq. + 1 
1 


{dV)Q. 

dV 

dV 


AX.V + -Cn 


dV 


where we have used Xq = 1 ■ Einally, we may write 


ds^ = —du^ V ^ du — Tr 
4 2 


dV- 


1 


A£,v + -Cn 


-dV 


1. Singular case 

In the singular case, i.e. when some of the eigenvalues of 
Vy are +1, the sum in dZTb is performed only along the ele¬ 
ments where WiWj + 1. The proof of this fact closely fol¬ 
lows an analogous observation in the fermionic case J^. 
Let ly = Zii w,|/){/| be the eigenvalue decomposition of W, 
where |i) is the eigenvector of ly with eigenvalue w; and 
let Ci - w“' 6 [-1,1], Ci - tanh(g, 72 ), where g,- are the 
symplectic eigenvalues of G. Using this notation, dW = 
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2/(1 - wj)^\i){i\ + Wi(\i){di\ + Inserting the above 

expression in (l27l l we find 


d Tr 


dW- 


-dW 


-Cw - 1 

Zt 


(Wi-Wj)^ T 

^\{di\j}\^ 


i*i 


The first term in the above equation is well-defined also when 
w/ —> +1. To prove that the second term is bounded we define 
f{x,y) = (x - y)^(l - xyy^ and write 


<5 = ^ ^(1 - wf)dgj -I- Yjf(ci,Cj)wiWj\{di\j)f- . (E3) 

i <*i 

As shown in Lemma 3 of Refi ll^ , the function f(x,y) is 
bounded in [-1,1]^ f{x,y) < 4, and lim(;,_ 3 ,)^(±i_±i)/(x,y) = 
0. Therefore, the elements such that w/Wy = 1 do not con¬ 


tribute in the sum (IE31 l. Numerically, this corresponds to tak¬ 
ing the pseudo-inverse of the superoperator in (IZTT i or, equiv¬ 
alently, in manually avoiding the sum over the elements such 
thatw/Wj = 1. Therefore, even though Ea. (IE3b has been found 
assuming that w, ^ +1, it can be analytically extended to the 
general case. 

Notice that for pure states, where w, = +1, the effect of the 
function / in (IE3b can be obtained equivalently by the function 
f(x,y) - {x — y)^/2. Taking this substitution in (IE3b we find 


dpure = '^f{Ci,Cj)WiWj\{di\i)\^ = ^ 


{Wi-WjY^, .,,2 


i*j 

= -Tr 
2 


IWjWj 


-\{di\j)r 


1 1 
— dW — dW 
IT IT 


(E4) 


The above equation provides a simpler expression for the Bu¬ 
res metric for a pure Gaussian state. 









